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INTRODUCTION 


Active  sonar  surveillance  systems  typically  experience  diminished  detection  performance  in 
reverberation-limited  environments.  In  some  underwater  environments,  the  transmission  of  an 
active  sonar  pulse  leads  to  acoustic  reflections  from  the  ocean  bottom,  from  the  sea  surface,  and 
from  within  the  ocean  volume.  These  reflections  can  have  significant  power  levels  at  the  receiver 
and  may  present  a  background  level  that  is  higher  than  the  level  of  the  local  ambient  noise.  This 
is  known  as  a  reverberation-limited  environment,  because  it  is  the  reverberation  and  not  the 
ambient  noise  that  limits  detection  performancel^l.  In  such  an  environment,  target  echoes  may 
not  be  strong  enough  to  be  detectable  above  the  reverberation  level,  being  masked  by  it.  The 
amount  of  reverberation  experienced  is  variable,  depending  largely  on  the  oceanographic  condi¬ 
tions,  acoustic  propagation,  and  surrounding  bathymetry.  In  addition,  reverberation  levels  are 
dependent  on  source  levels  and  the  transmitted  pulse  characteristics,  such  as  center  frequency, 
bandwidth,  and  duration.  There  are  a  variety  of  approaches  and  methodologies  that  reduce  the 
adverse  effect  of  reverberation,  enabling  the  detection  of  targets  of  interest. 

In  an  active  sonar  system,  the  effect  of  reverberation  can  be  reduced  by  either  the  spatial  or 
temporal  processing  of  the  received  signals.  In  spatial  processing,  adaptive  beamforming  can 
reduce  signal  leakage  through  the  beams’  sidelobes  that  are  due  to  strong  reverberation  features. 
This  is  done  by  steering  nulls  in  the  directions  of  the  interferers.  Reverberation  reduction  by  tem¬ 
poral  processing  is  based  on  the  ability  to  discriminate  a  target  signature  from  the  low  Doppler 
reverberation  field.  This  may  be  done  with  the  use  of  Doppler-sensitive  waveforms,  adaptive 
Doppler  filtering,  or  coherent  multiping  processing. 

The  temporal  multiping  processing  methods  used  to  reduce  reverberation  are  referred  to  as 
ping-to-ping  or  coherent  interping  reverberation-suppression  techniques.  These  methods  are 
based  on  the  assumption  that  the  characteristics  of  the  contaminating  reverberation  are  similar 
from  one  received  transmission  to  the  next.  If  this  assumption  holds,  then  by  using  other 
received  transmissions  as  reference  signals,  it  is  possible  to  subtract  out  the  common  reverbera¬ 
tion  they  share,  from  the  received  input  signal.  The  reference  signals  are  usually  taken  from  prior 
or  subsequently  received  pings.  The  subtraction  is  actually  implemented  with  an  adaptive  digital 
filter  in  order  to  optimize  the  reduction  of  the  reverberation. 

The  adaptive  suppression  algorithms  can  be  implemented  and  applied  to  datasets,  and  the 
resulting  performance  gains  can  be  measured.  Many  of  the  adaptive  algorithms  require  input 
parameters,  and  their  performance  is  sensitive  to  them.  In  addition,  the  performance  of  these 
techniques  depends  on  the  characteristics  of  the  acoustic  dataset  itself.  Therefore,  it  is  important 
to  be  able  to  determine  the  limits  of  performance  by  using  these  techniques  on  a  particular  data¬ 
set.  This  can  be  done  in  large  part  with  the  interping  coherence  measurement. 

The  interping  coherence  measurement  can  be  made  on  the  dataset,  and  from  it,  theoretical 
performance  gains  can  be  determined.  The  ability  to  predict  the  amount  of  gain  obtainable 
against  reverberation  is  valuable,  because  it  enables  the  assessment  of  the  effectiveness  of  these 
techniques  on  the  dataset.  Then,  if  the  data  are  processed  using  one  of  the  suppression 
algorithms,  the  parameters  may  be  adjusted  until  the  theoretical  gains  are  achieved. 
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While  it  is  possible  to  suppress  reverberation  using  these  techniques,  it  is  also  possible  that 
the  target  echo  undergoes  suppression  as  well.  For  these  techniques  to  be  effective,  it  is  essential 
that  the  target  echo  is  preserved,  while  only  the  reverberation  is  suppressed.  In  order  for  the  tar¬ 
get  echo  to  be  preserved,  it  must  decorrelate  from  ping  to  ping.  This  implies  some  target  motion. 
Therefore,  target  motion  relative  to  the  source  and  receiver  affects  the  suppression  techniques’ 
ability  to  preserve  target  echoes. 

This  report  addresses  the  performance  criteria  of  ping-to-ping  reverberation-suppression 
techniques.  First,  typical  sonar  signal  processing  methods  and  the  sonar  system  geometry  are 
reviewed.  Then,  the  interping  coherence  measurement  and  the  effects  of  target  motion  are 
addressed.  Finally,  interping  coherence  measurements  and  suppression  results  are  shown  for 
some  reverberation  data  collected  at  sea,  as  validation  of  the  algorithms’  performance. 

ACTIVE  SONAR  SOURCE 

Consider  an  acoustic  source  that  transmits  a  sequence  of  N  identical  signals,  as  shown  in 
figure  1.  Assume  that  each  individual  pulse,  p(r),  is  time  varying,  of  finite  duration,  Tp,  and  of 
bandwidth  B.  Let  the  ping  repetition  interval  be  Tprii  such  that  the  total  source  transmission 
sequence  is  given  by 

N  1  ^ 

Pse4(l)  =  X  H'  ■  M  (1) 

k=\  *=1 

where  >  Tp  and  pk(t)  is  the  ping  of  the  sequence.  The  transmitted  pings  are  assumed  to 
be  frequency  modulated  or  broadband  waveforms,  and  insensitive  to  Doppler  shifts.  These  types 
of  signals  are  used  because  of  their  excellent  range-resolution  properties.  Generally,  the  source 
array  has  vertical  aperture  to  provide  directionality,  and  it  transmits  in  all  azimuthal  directions. 
The  acoustic  transmissions  are  then  propagated  through  the  acoustic  channel.  The  transmissions 
may  experience  amplitude  fluctuations,  phase  shifts,  boundary  interactions,  and  propagation  over 
multiple  paths  by  the  time  they  arrive  at  the  receive  array. 


Figure  1.  Acoustic  ping  transmitting  sequence. 


2 


ACTIVE  SONAR  SIGNAL  PROCESSING 


Active  sonar  systems  typically  perform  signal  processing  functions  on  the  received  data  simi 
lar  to  those  shown  in  figure  2.  Each  block  in  the  processing  chain  is  summarized  below. 


DATACXXIECT'IN 

1  -SIGNAL 

TRANSOUCnOi 

-^-SIGNAL 
^  AMPUFtCATION 
-  SIGNAL  RLI^ING 
l-SIGNAL 

DIGmZATION 


COMPLEX 

OEMQOULATKM 

-  BANDPASS  FILTER 

-  SIGNAL 
BASEBANONQ 

-  DOWNSAMPLING 

-  COMPLEX  DATA 


BEAMFORMINQ 


-  SNR  GAIN 

-  SPATIAL  RLTER 


NORMAUZATION 

-  OBTAIN  SNR 

-  THRESHOLOING 

-  DATA  REDUCTION 

-  DETECTION 


REVERBERATION 

SUPPRESSION 


I  MATCHED  RLTER 

-  SNRQAIN 

-  PULSE 
COMPRESSION 


Figure  2.  Typical  active  sonar  signal  processing  path. 

DATA  COLLECTION 

Sonar  receive  arrays  are  composed  of  multiple  acoustic  sensors,  put  together  in  a  variety  of 
configurations,  spacings,  and  apertures.  Each  hydrophone  transducer  senses  the  pressure  incident 
upon  it  and  converts  the  acoustic  signal  into  a  voltage.  The  voltage  signals  are  then  passed 
through  the  receiver  electronics,  where  they  may  be  conditioned,  amplified,  and  filtered.  Follow¬ 
ing  this,  the  signals  are  digitized  with  the  use  of  an  analog-to-digital  converter.  With  the  data  in 
digital  form,  the  remaining  signal  processing  functions  can  be  efficiently  implemented  using 
digital  hardware  and  software. 

COMPLEX  DEMODULATION 

In  active  sonar,  except  for  small  Doppler  shifts,  returning  echoes  occupy  roughly  the  same 
frequency  band  as  the  transmitted  ping.  Since  everything  of  interest  is  within  this  band,  it  is 
desirable  to  filter  around  it.  Bandpass  filtering  over  the  frequency  extent  of  the  waveform 
effectively  eliminates  the  contamination  of  out-of-band  noise,  thereby  increasing  the  signal-to- 
noise  ratio.  Once  the  signal  frequency  band  is  passed,  it  is  possible  to  demodulate  the  data,  from 
the  ping’s  center  frequency  down  to  zero.  The  resulting  basebanded  data  are  complex-valued, 
which  is  a  convenient  form  for  the  signal  processing  fimctions  that  follow.  In  addition, 
demodulation  makes  it  possible  to  reduce  the  input  data  sample  rate  down  to  the  new  Nyquist 
frequency,  which  corresponds  to  the  bandwidth  passed  by  the  filter.  This  filtering  is  similarly 
done  on  all  the  available  hydrophone  element  data. 
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BEAMFORMING 


With  an  array  of  sensors  it  is  advantageous  to  sum  all  the  available  hydrophone  signals,  with 
appropriate  delays,  thereby  enhancing  the  signal-to-noise  ratio.  This  beamforming  process  also 
provides  the  spatial  discrimination  of  a  signal’s  arrival,  which  is  unobtainable  with  a  single 
omnidirectional  hydrophone.  Multiple  beams  can  be  formed  simultaneously,  and  steered  to  look 
in  specified  directions,  to  provide  the  desired  spatial  coverage. 

MATCHED  FILTER 

Because  the  transmitted  pulse  is  known,  active  sonar  processors  typically  make  use  of  the 
matched  filter.  A  matched  filter  correlates  the  received  data  stream  with  a  replica  of  the  trans¬ 
mitted  pulse.  The  matched  filter  provides  a  signal-to-noise  ratio  gain  when  the  replica  matches 
up  with  similar  signals  in  the  received  data  stream.  In  addition,  when  frequency-modulated 
broadband  pulses  are  used,  the  matched  filter  compresses  any  echoes  into  sharp  peaks.  This 
pulse-compression  property  is  desirable,  because  it  allows  for  the  resolution  of  echoes  at  time 
differences  equal  to  the  reciprocal  of  the  bandwidth.  In  addition,  any  effective  Doppler  shift  in 
the  echo  will  translate  directly  into  a  time  shift  of  the  correlation  peak.  The  matched  filtering  is 
performed  on  each  of  the  output  beams.  At  the  output  of  the  matched  filter,  target  echoes  may  be 
searched  for  on  the  square  of  the  magnitude  of  the  complex  data,  which  is  the  signal’s  power. 

REVERBERATION  SUPPRESSION 

At  this  stage  of  the  processing,  although  the  signal-to-noise  ratio  has  been  significantly 
increased,  target  echoes  may  or  may  not  be  detectable  above  the  noise  or  reverberation  back¬ 
ground  level.  If  the  desired  target  echo  is  above  the  ambient  noise  background,  but  obscured  by  a 
distributed  reverberation  background,  then  further  processing  may  improve  detection  perfor¬ 
mance.  In  such  an  event,  interping  reverberation  subtraction  techniques  may  be  employed  to  an 
advantage,  by  suppressing  the  obscuring  reverberation.  Multiple  pings  are  collected  and  then 
input  into  an  adaptive  interference-cancellation  filter. 

NORMALIZATION,  THRESHOLDING,  AND  DETECTION 

The  processed  data  are  normalized  as  a  part  of  detection  processing.  The  normalizer  outputs 
a  sequence  from  which  the  signal  power  to  local  background  (noise  or  reverberation)  power  ratio 
can  be  obtained.  This  ratio  is  a  measure  of  the  echo’s  detectability.  In  the  noise-limited  case, 
these  values  can  be  thought  of  as  the  signal-to-noise  ratio.  A  threshold  is  then  applied  to  the 
normalized  data,  with  only  the  data  that  are  above  the  threshold  passed  on  to  the  detector.  The 
detector  may  either  be  an  automated  computer  algorithm  providing  alerts  or  a  human  operator 
using  graphical  displays.  Detection  processing  is  performed  on  each  of  the  matched  filter  beam 
outputs. 


SONAR  SYSTEM  GEOMETRY 

Consider  the  bistatic  sonar  geometry  as  described  by  Coxl^l,  and  as  shown  in  figure  3.  A 
bistatic  geometry  is  the  most  general  one,  in  which  the  active  source  and  receiver  are  in  different 
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TARGET 


SOURCE  d„  RECEIVER 


Figure  3.  Sonar  system  geometry  (bistatic). 

physical  locations.  The  transmitted  signals  travel  from  the  source  to  the  receiver  along  the  direct 
path  dsr-  The  reception  of  this  signal  in  the  received  data  stream  is  called  the  “direct  blast,”  and  it 
typically  has  very  high  power,  since  the  acoustic  energy  is  transmitted  with  very  high  power  and 
undergoes  only  one-way  attenuation.  Assume  that  the  transmitted  pulse  also  travels  to  a  target  of 
interest  along  the  path  is  reflected  from  the  target,  and  travels  along  the  path  4>-  to  the 
receiver.  The  total  target  range  in  this  configuration  is  dtot,  which  is  equal  to  the  sum  of  <4/  and 
d,r.  The  propagation  path  for  the  target  echo  is  longer  than  the  direct  path,  by  an  additional  dis¬ 
tance  de,  given  as 

dg  —  df0i  ~  dgf  —  d^i  +  dfj,  —  d^  (2) 

These  ranges  can  also  be  described  equivalently,  in  terms  of  their  acoustic  travel  times,  as 

~  ^sl  —  tsr  (3) 

where  t  —  ~  and  c  is  the  velocity  of  sound  propagation  in  the  ocean.  The  time  tg  is  the  time  after 
the  arrival  of  the  direct  blast  that  a  target  echo  arrives  at  the  receiver.  The  source  and  the  receiver 
positions  form  the  foci  of  “equitime”  ellipses,  which  are  loci  of  points  for  which  the  bistatic 
travel  times  are  equal.  The  size  of  the  ellipse  corresponds  directly  to  the  echo  arrival  time,  tg. 

The  equitime  ellipse  for  a  source-receiver  pair  is  shown  in  figure  4.  In  order  for  the  echo  to  be 
cleanly  received  without  being  masked  by  the  direct  blast  energy  requires  that  Tp<tg.  Echoes 
from  scatterers  located  inside  the  equitime  ellipse  that  corresponds  to  Tp  are  drowned  out  by  the 
enormous  energy  of  the  direct  blast.  This  area  is  called  the  direct  blast  blanking  region. 
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REVERBERATOR 


TARGET 


ELUPSE 

Figure  4.  Equitime  ellipses  and  blanking  region. 

Typically,  the  bistatic  receiver  is  a  horizontal  line  array  of  hydrophones,  either  fixed  on  the 
bottom  or  towed  by  a  moving  platform.  It  typically  has  a  beamforming  capability  to  spatially 
filter  the  received  data.  Beamforming  enhances  the  ability  to  determine  the  arrival  angle  of 
signals  and  yields  a  gain  in  signal-to-noise  ratio.  Figure  4  also  shows  one  such  beam,  pointing  in 
the  direction  of  a  target  situated  in  a  highly  reverberant  area  populated  by  a  variety  of  scatterers. 
In  general,  for  a  linear  array,  there  is  an  ambiguity  to  the  beam  steering  directions,  and  the  inter¬ 
fering  reverberation  could  be  coming  in  from  the  front  or  the  back  beam  direction.  In  addition, 
strong  reverberation  from  other  directions  may  leak  through  the  beam’s  sidelobes,  degrading 
detection  performance. 

The  beam  directions  and  time  delays  directly  map  to  azimuthal  directions  and  equitime 
ellipses,  respectively.  From  the  processed  data  and  this  mapping,  an  estimate  of  the  approximate 
physical  location  of  the  echo-producing  target  can  be  made.  Figure  5  depicts  the  magnitude  of  a 
typical  matched  filter  power  output  for  a  target-pointing  beam,  for  a  noise-limited  case.  Note  the 
reception  of  the  direct  blast,  and  at  time  te  later,  the  reception  of  the  target  echo. 

The  sonar  geometry  has  been  described  for  the  bistatic  configuration.  This  case  is  empha¬ 
sized  here,  since  it  is  the  most  general  case  for  a  sonar  system.  The  more  familiar  configuration 
is  monostatic,  where  the  source  and  the  receiver  are  colocated.  The  monostatic  geometry  is  just  a 
special  case  of  the  bistatic  configuration.  This  occu''s  when  d^r  =  0,  and  dsi  =  dtr  =  4-  Therefore, 
for  the  monostatic  case,  ^4  =  24-  Also,  the  equitime  ellipses  lose  their  eccentricity  and  become 
circles. 
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Figure  5.  Typical  matched  filter  output,  noise-limited, 

REVERBERATION-LIMITED  ENVIRONMENT 

To  ensure  that  received  reverberation  is  similar  from  ping  to  ping,  only  the  case  where  both 
the  source  and  the  receiver  positions  remain  stationary  is  considered.  This  kind  of  bistatic  sonar 
system  has  a  “fixed-fixed”  configuration.  The  fixed-fixed  configuration  increases  the  likelihood 
of  having  stationary  and  coherent  reverberation  in  the  received  data  stream  from  one  transmis¬ 
sion  to  another.  The  receive  array  receives  the  direct-blast  arrival  from  the  transmitted  pulse 
Pk{t),  as  well  as  any  scattered  reflections  from  the  bottom,  from  the  surface,  and  from  within  the 
ocean  volume.  In  this  fixed-fixed  configuration,  the  reverberation  echoes  associated  with  each 
transmission  have  the  best  chance  of  being  coherent  from  ping  to  ping.  Coherence  between  the 
reverberation  from  one  ping  to  the  next  is  a  prerequisite  for  the  subtraction  algorithms  to  be 
effective,  and  will  be  demonstrated  later  in  this  report. 

The  received  signal  from  the  it:**’  transmitted  pulse  pk,  can  be  characterized  as  the  sum  of 
all  the  individual  reflections  from  each  scatterer,  plus  the  local  ambient  noise  background.  Each 
reverberation  feature  will  have  its  own  scattering  strength  and  will  yield  a  distributed  collection 
of  echoes  in  the  received  data  stream  at  the  appropriate  times  and  in  the  appropriate  beams. 
Assuming  a  stationary  reverberation  field  with  identical  (perfectly  coherent)  reverberation  from 
ping  to  ping,  each  received  signal  Xk  can  be  given  as 

+  «*(0  (4) 

where  is  the  same  coherent  reverberation  component  and  n^  is  the  differing  local  noise  signal 
for  each  ping. 
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Now,  consider  the  received  reverberation  power  caused  by  scattering  from  the  ocean  bottom. 
Pronounced  bathymetric  features,  continental  shelves,  and  the  ocean  bottom  in  shallow  water 
tend  to  be  very  good  scatterers  of  sound,  resulting  in  an  overall  increase  in  the  received  signal 
power  above  that  of  the  ambient  background  noise  power.  It  is  this  type  of  reverberation  to 
which  interping  reverberation-suppression  techniques  are  best  suited,  because  the  properties  and 
statistics  of  the  reverberation  are  most  likely  to  be  the  same  from  transmission  to  transmission. 

The  increases  in  received  signal  power  caused  by  bottom  reverberation  occur  in  the  time 
series  data  in  the  same  manner  as  explained  previously  for  target  ecb'tes.  The  time  cell  of  the 
received  reverberation  corresponds  to  an  equitime  ellijsse.  The  azimuthal  direction  of  the  scat¬ 
tered  energy  may  be  determined  by  the  formation  of  multiple  beams.  A  depiction  of  a  single 
beam’s  matched  filter  power  output  with  reverberation  features  might  appear  as  in  figure  6, 
where  the  increases  in  signal  power  are  due  to  bathymetric  features.  Notice  that  when  the  rever¬ 
beration  power  rises  above  the  ambient  noise  power,  the  sonar  system  is  reverberation-limited  in 
those  regions.  In  this  situation  target  echoes  may  be  masked  and  therefore  remain  undetectable. 
With  these  conditions,  the  suppression  techniques  may  provide  detection  performance  gains  by 
reducing  the  amount  of  obscuring  reverberation. 

COHERENT  REVERBERATION  SUPPRESSION 

Consider  a  target  that  is  moving  with  velocity  vj  within  an  area  of  high  bottom  reverberation 
and  obscured  by  it  in  the  matched  filter  output,  as  shown  in  figure  6.  In  the  absence  of  reverbera¬ 
tion,  the  target  echo  is  detected  easily,  since  it  is  well  above  the  ambient  noise  background.  Here, 
however,  the  reverberation  power  obscures  it.  Assume  that  a  ping  sequence  PseqiO  is  transmitted 
with  the  same  waveform  type  and  with  ping  repetition  interval  Tpri-  The  repetition  interval  is 
typically  chosen  so  that  the  reverberation  has  diminished  significantly,  to  levels  below  that  of  the 
ambient  noise  background,  before  the  next  ping  is  received.  This  eliminates  undesirable  rever¬ 
beration  buildup  from  one  ping  to  the  next.  In  addition,  there  may  be  other  factors  driving  this 
choice,  such  as  source  duty  cycle  limitations.  The  acoustic  energy  propagates  to  the  target  and 
other  scatterers,  reflects  from  them,  and  propagates  to  the  receiver.  Figures  7a  and  7b  depict  two 
consecutively  received  signals  at  the  output  of  a  matched  filter.  These  pings  have  been  extracted 
from  the  received  data  stream  and  time-aligned  by  Tpri  so  that  the  arrival  of  their  direct  blasts 
line  up.  Comparing  the  two  pings,  it  is  seen  that  the  reverberation  power  appears  to  be  similar 
from  ping-to-ping.  The  reverberation  processing  is  essentially  an  adaptive  subtraction  of  one 
ping  from  another.  If  the  reverberation  background  is  coherent,  then  the  suppression  techniques 
may  be  effectively  applied.  Because  of  the  target’s  motion,  its  echoes  shift  in  arrival  time  relative 
to  one  another  from  one  ping  to  the  next.  With  sufficient  time  shift  between  the  echoes,  the  input 
target  echo  can  be  passed  through  the  reverberation  processing  unsuppressed.  An  example  output 
of  the  suppression  process  is  depicted  in  figure  7c.  Notice  that  with  the  reduction  in  reverbera¬ 
tion,  the  target  has  become  detectable. 
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NdSE-UMITEO 

REVERBERATION-UMITED 


Each  individual  received  ping’s  matched  filter  output  can  be  thought  of  as  a  single  record  in  a 
collection  or  ensemble  of  records  forming  the  random  variable  Xi{t).  Each  ping  is  aligned  in  the 
ensemble  so  that  the  beginning  of  each  subsequent  record  is  some  multiple  of  after  the 
beginning  of  the  current  record.  Done  this  way,  the  direct  blasts  and  reverberation  features  line 
up.  Xi(t)  is  not  ergodic  or  stationary,  in  the  strictest  sense,  because  the  ensemble’s  mean  received 
power  level  may  vary  with  time.  This  is  because  changing  reverberation  levels  are  received 
through  time  from  reverberation  features  at  different  bistatic  ranges.  However,  the  statistics  for  a 
short  localized  time  window  bt,  centered  at  an  arbitrary  time  t',  may  be  similar  between  two  or 
more  records  of  the  ensemble.  For  example,  the  power  level  may  vary  similarly  through  time  for 
one  or  more  records  of  Xf.  This  is  also  shown  in  figure  7a,b.  Assuming  similar  acoustic  propaga¬ 
tion  conditions,  constant  source  level,  and  a  Exed-fixed  geometry,  localized  reverberation 
regions  of  the  time  series  should  have  similar  power  levels.  Likewise,  higher  order  statistics  may 
also  be,  in  this  sense,  “stationary”  from  record  to  record  across  the  ensemble  at  a  particular  time 
location  t'.  The  interping  coherence  is  a  measurement  that  is  made  between  pairs  of  records  in 
the  ensemble,  over  the  time  interval  dt.  This  measurement  is  made  repeatedly  at  periodic 
intervals  along  the  time  axis,  resulting  in  the  interping  coherence  function.  An  example  of  the 
interping  coherence  function  is  shown  in  figure  7d.  From  values  of  interping  coherence,  the 
theoretical  amount  of  reverberation  reduction  can  be  determined,  as  will  be  described  in  a  later 
section. 
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Figure  7.  Reverberation  processing  and  interping  coherence. 
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ADAPTIVE  REVERBERATION-SUPPRESSION  ALGORITHMS 


The  ping-to-ping  subtraction  description  given  is  somewhat  simplified.  Straight  coherent 
subtraction  of  one  ping  from  another  rarely  yields  the  desired  output,  as  shown  in  figure  7c.  This 
is  because  the  two  received  signals  are  not  identical,  since  there  may  be  amplitude  variations, 
time  delays,  and  phase  shifts  between  the  reverberation  in  each  of  them,  as  well  as  a  differing 
local  noise  component.  The  actual  suppression  algorithms  are  implemented  with  an  adaptive 
filter  structure,  so  that  an  optimal  estimate  is  obtained  for  the  reference.  Subtraction  of  this 
estimate  from  the  input  ping  yields  an  output  time  series  with  minimum  power.  The  degree  to 
which  the  adaptive  algorithms  can  do  this  is  determined  by  the  coherence  between  the  two  (or 
more)  time  series  going  into  the  process. 

The  principles  of  using  an  adaptive  filter  structure  for  the  cancellation  of  interference,  along 
with  some  applications  in  the  areas  of  electrocardiography,  speech,  and  anays,  have  been 
described  in  detail  by  Widrow  et  al.f^l.  If  a  correlated  reference  of  the  interference  is  available,  it 
can  be  used  to  obtain  an  optimum  subtraction  of  the  interference  in  the  input.  Figure  8  shows  a 
block  diagram  representing  the  adaptive  filter  configuration.  The  filtered  reference  signal  is  sub¬ 
tracted  from  the  input  to  obtain  the  output.  The  goal  is  to  minimize  the  power  of  the  output  sig¬ 
nal.  To  do  this,  the  filter  weights  are  adaptively  updated  as  the  data  are  passed  through  it,  so  that 
the  output  yields  the  least  possible  power.  The  adaptive  weight  control  mechanism  may  be  one  of 
several  types  (e.g.,  least  mean  square  (LMS)  and  recursive  least  squares  (RLS));  generally  it  uses 
prior  values  of  the  output  to  modify  the  current  set  of  filter  weights.  In  this  manner,  the  feedback 
loop  drives  the  filter  to  track  the  variations  in  the  input  data,  as  well  as  keep  the  output  power  of 
the  filter  as  small  as  possible.  Figure  8  shows  that  the  reference  may  be  composed  as  the  output 
of  a  parallel  bank  of  adaptive  filters,  each  operating  on  a  single  reference  channel.  In  this 
application,  the  input  and  the  references  are  obtained  from  sequentially  received  pings.  Notice 
that  in  this  configuration,  the  input  ping  is  not  included  in  the  reference,  which  is  referred  to  as 
Case  1. 


INPUT  (PING  1) 


REFERENCE  (PING  2) 
REFERENCE  (PING  3) 


■ 


■ 


REFERENCE  (PtNG  N) 


Figure  8.  Adaptive  suppression  using  LMS  or  RLS  (Case  1). 
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There  are  a  variety  of  algorithms  used  to  adaptively  update  the  filter  weights  for  optimum 
performance.  Detailed  descriptions  of  these  algorithms,  and  their  performance  characteristics, 
are  given  by  Haykin,  along  with  an  extensive  reference  listl'^l  Two  of  the  more  common  algo¬ 
rithms  are  the  LMS  and  RLS  algorithms.  Although  their  goal  is  the  same,  the  various  adaptation 
schemes  perform  differently,  and  one  may  be  better  suited  to  a  particular  application  than  others. 
The  differences  in  the  algorithms  usually  have  to  do  with  their  tracking  performance,  their 
ability  to  adapt  and  converge  to  the  optimum  filtering  solution,  and  the  computational  require¬ 
ments  and  complexity  of  the  filter  implementation. 

The  adaptive  algorithms  generally  have  at  least  two  parameters  that  must  be  specified:  the 
adaptation  rate  or  memory  and  the  number  of  taps  in  the  filter.  Although  it  is  not  our  intent  here 
to  discuss  the  details  of  these  parameters,  it  is  important  that  appropriate  choices  of  parameters 
be  made,  choices  that  are  well  suited  or  matched  to  the  data  being  processed.  One  parameter 
controls  the  rate  at  which  the  adaptive  filter  will  update,  in  its  convergence  to  the  optimal  filter 
solution.  It  is  necessary  that  the  filter  adapt  fast  enough  to  track  the  changes  in  the  distributed 
reverberation  to  be  suppressed,  but  not  so  fast  that  it  also  removes  the  target  echo.  In  addition, 
the  number  of  adaptive  weights  should  be  selected  to  achieve  optimum  performance  for  the  data 
being  processed. 

Figure  9  shows  an  alternative  approach  for  multiple  ping  reverberation  suppression.  Here,  the 
sequentially  received  pings  are  considered  to  be  a  rectangular  data  matrix,  with  the  input  ping 
also  included.  A  singular  value  decomposition  (SVD)  is  then  performed  on  the  data  matrix.  The 
SVD  is  essentially  an  eigenvalue,  eigenvector  decomposition  of  the  data  matrix.  The  output  of 
the  SVD  is  an  estimate  of  the  input,  reconstructed  with  different  eigenvalue-weighted  combina¬ 
tions  of  the  eigenvectors.  Assuming  that  the  reverberation  is  coherent  from  ping  to  ping,  the 
strongest  eigenvector  is  composed  essentially  of  the  reverberation,  the  principal  component  to  be 
removed.  Therefore,  the  output  of  the  suppression  processing  is  the  input  ping  with  one  or  more 
of  the  strongest  eigenvectors  removed.  Figure  9  shows  the  case  where  the  output  is  the  input 
ping  with  the  single  strongest  eigenvector  removed.  This  technique  does  not  require  any  input 
parameters,  but  at  the  output  there  may  be  some  uncertainty  as  to  which  eigenvector  or  eigen¬ 
vectors  to  remove.  This  approach  is  considered  as  Case  2,  differing  from  Case  1  because  it 
includes  the  input  ping  in  the  makeup  of  the  reference.  It  will  be  shown  that  while  this  technique 
reduces  more  reverberation  than  the  Case  1  technique,  it  suffers  an  additional  loss,  so  that  its 
overall  performance  against  reverberation  is  equivalent  to  that  of  the  Case  1  technique. 

Because  the  performance  of  the  adaptive  algorithms  is  sensitive  to  the  accompanying 
parameters  and  the  particular  dataset  the  algorithm  are  of>erating  on,  a  reliable  method  to  assess 
the  theoretical  performance  gains  is  desired.  Knowing  the  bounds  of  performance  we  can  then 
tune  the  algorithms’  parameters  until  that  bound  is  approached.  The  theoretical  performance 
criteria  are  now  discussed. 
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Figure  9.  Adaptive  suppression  using  SVD  (Case  2). 


INTERPING  COHERENCE 


The  normalized  cross-correlation  function  p,  between  two  real-valued  signals  ac^(r)  and 
is  givenl^l  by  the  cross-correlation  function  divided  by  the  power  in  the  two  input  signals,  as 


SB 

I 


Xw(t)  •  x„(t  -  t)dt 


00  00  " 

I  j  MOP* 


1/2 


(5) 


With  this  normalized  form  for  correlation,  p  can  take  on  values  between  -1  and  +1.  Now,  con¬ 
sider  the  case  of  the  two  signals  and  j(^,  received  from  two  different  transmissions,/!^  andp^. 
The  signals  are  so  aligned  by  their  direct-blast  arrival  that  the  reverberation  features  line  up,  and 
a  short  time  window  of  data,  6t,  is  obtained  from  each  received  signal.  Using  a  discrete  time 
implementation,  the  normalized  cross-correlation  function  is  estimated  by  using  the  two  Hnite- 
length  segments.  Assume  that  the  maximum  value  of  the  estimate  of  p„^„  occurs  at  t  =  ‘^,  the 
delay  lag  at  which  the  two  signals  are  most  correlated.  The  range  of  lag  values  for  x  is  from 
-Xmaxl'2-  to  +Tinfl3c/2,  where  Xmax  may  be  specified  so  that  bounds  are  placed  on  the  correlation 
search,  if  desired.  The  interping  coherence  estimate  ^m,n  is  defined  as  the  absolute  value  of  Pm^„ 
at  the  delay  x  =  1^,  such  that 

rni,n  ~  I  max[pin,n  (x)|  I  =  |pm,n(tp)|  (6) 


Here,  is  bounded  by  0  and  1  and  measures  the  similarity  or  coherence  between  the  two 
signals.  This  correlation  process  is  then  repeated  at  periodic  intervals  along  the  time  base  of  the 
received  pings.  Each  correlation  window  overlaps  the  previous  one  by  a  selectable  amount,  6o/. 
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This  results  in  the  interping  coherence  as  a  function  of  time,  r„,„(r),  which  is  also  dependent  on 
the  parameters  dr,  6o^,  and  Xmax-  interping  coherence  function  will  indicate  any  regions  in 
range  where  the  reverberation  is  similar  from  ping  to  ping. 

Now,  consider  the  interping  coherence  measurements  for  a  single  time  window,  ht,  in  terms 
of  the  reverberation  and  noise  that  make  up  each  signal.  For  finite-length  discrete  signals,  the 
correlation  integral  is  a  time-averaged  product,  which  is  denoted  by  the  symbol  <>.  In  addition, 
real-valued  signals  are  considered,  and  the  time  dependence  is  assumed  and  dropped  for 
simplicity  in  notation.  The  following  assumptions  are  made: 

1.  Each  received  ping,  and  Jt^,  is  comprised  of  a  reverberation  component  and  a  local 
ambient  noise  component  as 

“  ^»n  ■*’#1  ~  (7) 

2.  The  reverberation  components  in  the  and  pings  are  identical,  perfectly 
coherent,  and  given  by  r^.  In  addition,  the  coherent  power.  Pc,  is  given  by  the  averaged 
square  of  the  reverberation  component: 

rm  =  r„=  Tc 

Pc  =  (>•?)  (8) 

3.  The  coherent  reverberation  component,  r^,  is  uncorrelated  with  the  ambient  noise  in 
each  ping  so  that  their  averaged  product  is  zero: 

(rc  •  /«m)  =  {rc  •  n„)  =  0  (9) 

4.  The  noise  process  in  each  signal  is  random,  white  over  the  signal  frequency  band,  and 
uncorrelated  with  each  other.  In  addition,  the  averaged  noise  power,  P^,  is  assumed  to 
be  the  same  in  each  signal: 

(«m  •  «„)  =  0 

=  {nl)  =  =  Ps  (10) 

5.  The  interping  coherence  is  equivalent  between  all  pairs  of  pings  going  into  the 
reverberation-suppression  processing: 

Tm,/*  ~  r  (11) 

for  all  interping  combinations  m,  n. 


If  these  foregoing  assumptions  are  made,  the  interping  coherence  may  be  described  in  term 
of  the  reverberation  and  noise  power,  as  follows;  The  interping  coherence  between  the  /n*  and 
pings  will  then  be  given  as 


{ire  +  rim)  '  jrc  +  n„)) 
iirc  +  «m)^>  •  {(rc  +  n„)^) 


(r|  +  rciim  +  rcn„  +  n„n„) 
l{rl  +  2rcn„  +  n^)  •  (r|  -I-  2rcn„  +  nl) 


If  assumptions  2, 3,  and  4  are  made,  the  reveiberation  and  noise  cross-terms  are  zero,  so  that 
the  interping  coherence  simplifies  to 
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where  Pc  is  the  power  of  the  coherent  reverberation  component,  Pn  is  the  noise  power,  and  Pj  is 
the  total  power  of  both.  Therefore,  the  interping  coherence  is  equal  to  the  ratio  of  the  coherent 
reverberation  power  to  the  total  power  in  the  received  ping.  As  the  power  of  the  noise  increases, 
the  interping  coherence  decreases,  and  as  the  noise  decreases,  the  coherence  improves.  The 
interping  coherence  function  may  be  estimated  by  measurement,  as  in  Eq.  6.  This  measurement 
can  be  made  for  every  pair  of  pings  going  into  the  process.  An  example  of  the  interping  coher¬ 
ence  function  for  two  pings  is  depicted  in  figure  7d.  In  reality,  for  the  noise-only  regions,  the 
estimated  coherence  function  does  not  drop  to  zero  as  the  expression  for  coherence  predicts.  This 
is  because  there  may  be  some  small  correlation  due  to  the  finite  length  of  the  data  going  into  the 
estimation  of  the  conelation.  Because  of  this,  the  estimation  of  the  coherence  function  may 
exhibit  a  slightly  higher  floor  than  that  for  a  completely  uncorrelated  signal.  In  addition,  the 
cross-correlation  function  has  freedom  to  search  over  some  range  in  x  for  possible  correlations. 
Hence,  although  the  noise  signals  are  uncorrelated,  there  may  be  some  significant  nonzero  cor¬ 
relation  at  other  delays,  which  is  picked  as  the  peak.  In  any  case,  the  levels  are  typically  still 
smaller  than  0.5  for  appropriate  choices  of  %„ax‘ 

Finally,  the  power  of  the  coherent  reverberation  component  can  be  expressed  as  a  function  of 
the  noise  power  and  the  interping  coherence  from  Eq.  13  as 

Pc  =  ^  (14) 

From  these  results,  expressions  can  be  obtained  for  the  theoretical  amount  of  reverberation  that 
may  be  reduced.  Two  cases  are  considered.  In  Case  1,  the  reference  ping  to  be  subtracted  is  made 
up  of  one  or  more  pings  other  than  the  input  ping.  In  Case  2,  the  input  ping  itself  is  included  in 
the  formulation  of  the  reference.  When  all  suppression  gains  and  echo  losses  are  taken  into 
account,  the  performances  of  Case  1  and  Case  2  implementations  are  equivalent,  when  the  same 
number  of  pings  are  used.  It  will  be  shown  that  while  Case  2  suppresses  more  reverberation  than 
Case  1,  it  also  suffers  an  additional  loss  in  target  echo  level,  so  that  its  overall  performance  is 
equivalent. 


REVERBERATION  POWER  REDUCTION  (CASE  1) 

Consider  the  reverberation  subtraction  process  as  depicted  in  figure  8.  A  total  of  V  pings  are 
to  be  used  in  the  reverberation-suppression  processing.  One  of  the  pings  is  chosen  as  the  input 
ping,  and  the  remaining  pings  are  used  in  the  formulation  of  the  reference  ping.  When  the  input 
ping  is  not  included,  the  reference  may  be  modeled  as  the  average  of  M  pings,  where  M  =  N-1. 
The  output  of  the  suppression  processing,  y,  is  the  reference  subtracted  from  the  input  as 

y  =  Xi„vut  -  \cf  (15) 

The  reference  ping  is  modeled  as  the  average  of  Af  pings  so  that 

M  M 
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Now,  from  Eq.  15  and  16,  the  output  is  given  by 
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=  n 
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i-l 


(17) 


input 


i«l 


Detection  processing  is  done  on  the  output  power  of  the  data.  The  output  power,  Py,  is 
obtained  by  squaring  the  output  signal  as 


P,  =  0,2)  =  „2  +  +  XJ 


(18) 


For  a  simple  two  ping  subtraction  case  (\f=l),  the  output  power  reduces  to  twice  the  noise 
power  of  the  input. 


Now,  from  Eq.  14  and  18,  the  reverberation  processing  power  gain  for  this  case,  ai,  can  be 
obtained.  The  gain  is  the  ratio  of  output  power  to  input  power  as 


Py  (i+AV’a'  (l+iiK 
P.  Pc  +  Pn  +  Pp 


_  (M  +  1)  ■  (1  -  D  ^  N(l-r)  ,J5,. 

M  (N-1) 

The  reciprocal  of  the  power  gain  gives  the  reverberation  power  reduction  (RPR),  defined  for 
Case  1  as 

RPR,  =  101og,„(i)  =  (20) 

Equation  20  gives  the  theoretical  amount  of  reverberation  power  (in  decibels)  that  is  reduc> 
ible,  as  a  function  of  the  interping  coherence.  As  the  measured  interping  coherence  increases,  the 
RPR  increases.  Figure  7c  depicts  the  power  output  of  a  two-ping  subtraction.  Notice  how  areas 
of  large  reduction  in  reverberation  correspond  to  areas  of  high  interping  coherence. 

Figure  10  shows  theoretical  plots  of  RPRi  as  a  function  of  F,  for  three  different  values  of  N. 
Notice  that  while  the  RPR  may  take  on  negative  values  for  small  values  of  interping  coherence, 
only  the  positive  portions  are  plotted.  This  is  the  applicable  region  of  the  curves,  since  the  adap¬ 
tive  algorithms  do  not  yield  output  power  greater  than  the  amount  in  the  input.  It  can  be  seen  that 
as  more  pings  are  used  in  the  reference,  greater  reductions  are  obtained  against  reverberation,  for 
a  fixed  coherence  value.  The  upper  bound  of  reduction  is  shown  for  =  «>,  which  is  a  constant  3 
dB  above  the  N  =  2  case,  for  any  given  coherence  value.  For  the  same  suppression  level,  the 
multiple-ping  reference  does  not  require  values  of  F  to  be  as  high,  particularly  at  smaller  RPR 
levels.  It  is  also  seen  that  large  reductions  in  reverberation  require  a  high  degree  of  interping 
coherence.  The  N  =2  case  requires  values  of  coherence  greater  than  0.5  to  achieve  any  rever¬ 
beration  reduction  at  all. 
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Figure  10.  Theoretical  RPR}  vs.  interping  coherence. 


REVERBERATION  POWER  REDUCTION  (CASE  2) 


In  Case  2,  the  input  signal  is  included  in  the  formulation  of  the  reference  signal.  This  is  the 
case  for  a  sequence  of  pings  input  into  an  SVD  algorithm,  as  shown  in  figure  9.  Again,  assume 
that  a  total  of  N  pings  are  used  in  the  reverberation  processing  and  that  one  of  the  pings  is  chosen 
as  the  input  ping.  Here,  the  reference  is  considered  to  be  an  average  of  the  total  of  all  N  pings, 
including  the  input,  such  that 


N 


Tr  + 


n 


input 


f  =  1 


N 


N-\ 


i=1 


(21) 


where  the  input  ping  is  taken  outside  of  the  summation.  As  before,  the  output  is  obtained  by  sub¬ 
tracting  the  reference  from  the  input  ping  as 


(22) 


The  output  power  of  the  suppression  processing  is  obtained  by  squaring  the  output  signal  as 


(23) 


Comparing  this  result  with  Case  1  (Eq.  18)  shows  a  reciprocal  relationship  between  the  factors 
(N  -  1)/N  and  (M  +  1)/M  multiplying  Pn-  This  is  because  in  this  case,  the  input  is  included  in  the 
reference,  whereas  in  Case  1,  it  is  not.  Consequently,  the  output  power  for  Case  2  will  always  be 
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smaller  than  for  Case  1  for  the  same  choice  of  N,  with  less  difference  as  ^  increases.  As  before, 
the  suppression  gain  is  the  ratio  of  output  power  to  input  power,  given  in  this  case  by 

Un-i)]p 

p,  [  y  f'<  (AT-D-d-n 
The  reverberation  power  reduction  for  Case  2  is  given  by 


N 


(24) 


(25) 


Plots  of  RPRi  are  shown  in  figure  11.  Here,  compared  with  Case  1,  higher  levels  of  rever¬ 
beration  reduction  are  possible.  While  this  is  indeed  true,  it  is  somewhat  misleading,  as  this 
method  also  yields  a  reduction  in  target  echo  level.  This  is  due  to  the  fact  that  the  target  echo  in 
the  input  ping  is  also  a  part  of  the  reference  and  therefore  partially  subtracts  itself  out.  This  is  the 
reference  loss  (RL),  which  only  applies  to  Case  2  implementations,  and  it  is  discussed  in  more 
detail  in  a  subsequent  section.  Notice  the  reciprocal  nature  of  the  curves  relative  to  Case  1.  The 
upper  bound  occurs  for  N=2,  and  the  lower  bound  is  when  N=<x>.  Also  notice  that  for  N=«), 
RFRi=RPR2. 


COHERENCE  OF  REVERBERATION 

The  interping  coherence  has  been  described  based  on  the  assumption  that  the  reverberation  is 
perfectly  coherent  between  the  pings  used  in  the  process.  In  reality,  this  may  not  be  so  in  many 
acoustic  environments.  In  addition  to  the  decrease  in  interping  coherence  due  to  uncorrelated 
noise,  there  may  also  be  some  decrease  due  to  reverberation  that  is  only  partially  coherent  This 
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effect,  along  with  the  impact  of  the  input  ping’s  reverberation-plus-noise-to-noise-power  ratio 
(RNNR),  is  now  explained.  The  RNNR  can  be  expressed  as 


RNNR  = 


N  _ 

Pn 


where  the  assumptions  made  for  interping  coherence  previously  still  apply.  This  quantity  can  be 
estimated  by  measuring  the  total  received  power  and  dividing  by  an  estimate  of  the  noise  power. 
The  noise  power  estimate  may  be  obtained  from  the  received  data  stream  just  before  active  trans¬ 
missions  or  by  measuring  the  noise  power  in  an  adjacent  frequency  band  just  outside  the  trans¬ 
mitted  signal.  The  RNNR  is  important  because  it  indicates  whether  the  data  are  noise-limited  or 
reverberation-limited.  For  a  reverberation-limited  case.  Pc  >  Pn,  and  the  RNNR  is  greater  than  3 
dB.  For  a  noise-limited  case.  Pc  <  Pn>  and  the  RNNR  is  less  than  3  dB.  When  no  reverberation 
is  present  (noise-only  case),  Pc  =  0,  and  the  RNNR  =  0  dB. 

Assuming  a  perfectly  coherent  reverberation  component,  the  maximum  possible  values  of  the 
interping  coherence  can  be  determined  by  measuring  or  estimating  the  RNNR.  Assume  that 
before  the  start  of  any  transmissions,  the  average  ambient  noise  power  level  is  Psambt  if 

remains  constant.  Once  the  transmissions  have  begun,  the  received  signal  power  is  the  sum  of 
the  reverberation  power.  Pc,  and  the  noise  power.  The  maximum  value  that  F  can  achieve,  based 
on  the  RNNR,  can  be  obtained  with  estimates  of  total  received  power  and  noise  power  as 

T'  _  Psamb  _  j  _  1 


Fmar  ~ 


RNNR 


The  maximum  interping  coherence  is  shown  in  figure  12,  as  a  function  of  the  input’s  RNNR.  It  is 
clear  from  this  figure  that  as  the  reverberation  level  drops  relative  to  the  noise,  the  maximum 
interping  coherence  drops,  resulting  in  less  suppression  performance.  But  because  there  is  less 
reverberation  in  the  input,  there  is  less  reverberation  power  to  suppress.  If  the  assumptions  out¬ 
lined  in  the  Interping  Coherence  section  are  indeed  true,  then  the  measured  interping  coherence 
will  be  equal  to  this  maximum  coherence  level.  In  that  case,  the  interping  coherence  is  only  a 
function  of  the  input  RNNR. 

Now,  consider  the  case  where  the  reverberation  component  is  not  perfectly  coherent  such  that 

Tn*  Tc  (28) 

Considering  only  the  reverberation  component,  in  the  absence  of  noise,  the  coherence  for  rever¬ 
beration  is 


reverb 


_  {rm  '  r„) 


The  interping  coherence  that  is  measured  using  cross-conelation  is  bounded  by  the  maxi¬ 
mum  coherence  and  reduced  by  the  coherence  of  the  reverberation  component,  as 

r  =  ^reverb  '  ^mwc  ( 
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Figure  12.  Maximum  coherence  obtainable  vs.  RNNR. 


Figure  13  shows  curves  of  the  RPRi  (N  -  2)  obtainable  as  a  function  of  input  RNNR  for 
various  values  of  Notice  that  at  high  RNNRs,  the  RPR  is  decreased  significantly  as  the 

coherence  of  the  reverberation  component  is  decreased  only  slightly.  Again,  as  before,  the  plot 
only  shows  the  curves  for  positive  RPRs,  since  the  algorithms  do  not  allow  output  power  to 
exceed  input  power.  Notice  that  in  this  case,  no  gains  can  be  made  if  the  RNNR  is  less  than  3  dB. 
This  means  that  only  when  Pc  >  Pn  ^re  any  gains  made  for  a  two-ping  suppression  process. 

Also,  the  uppermost  curve  (Tnverb  =  f  0)  shows  that  even  with  the  maximum  amount  of  coher¬ 
ence  in  the  reverberation  component,  there  is  still  a  residual  RNNR  of  3  dB  after  suppression. 
This  implies  that  no  reductions  of  reverberation  power.  Pc,  are  possible  below  the  power  of  the 
noise,  Ps,  for  a  two-ping  process.  The  coherence  of  the  reverberation  can  be  estimated  by  mea¬ 
suring  the  interping  coherence  function  and  dividing  by  the  estimate  of  Fmox  iti  Eq.  27.  F^ox  niay 
be  estimated  by  measuring  the  total  power  and  the  power  of  the  noise  as  described.  In  this 
manner,  the  coherence  of  the  received  reverberation  component  can  be  estimated. 

Figure  14  shows  a  similar  set  of  curves  for  Case  2,  with  five  pings  going  into  the  process. 
Here,  with  multiple  pings  used,  there  is  an  increase  in  reverberation  suppression  for  all  values  of 
Treverb-  Still,  as  before,  the  suppression  drops  as  the  coherence  of  reverberation  decreases.  For 
^reverb  =  the  RPR  is  1  dB  larger  than  the  input  RNNR.  Even  for  a  noise-only  case  (RNNR  =  0 
dB),  with  no  reverberation  present,  the  input  signal  is  reduced  by  1  dB.  This  is  because  the  input 
is  included  in  the  reference,  resulting  in  additional  suppression.  It  is  important  to  note  that  in  this 
case,  a  small  loss  in  both  noise  and  target  echo  levels  will  occur.  This  effect  is  described  in  detail 
in  a  later  section. 
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Figure  13.  Theoretical  RPR\  vs.  input  RNNR  (N  =  2). 


Figure  14.  Theoretical  RPR2  vs.  input  RNNR  (N  =  5). 

TARGET  ECHO  LOSSES 

While  the  ability  to  suppress  reverberation  is  indeed  valuable,  it  does  little  good  if  the  target 
echo  is  also  suppressed  or  eliminated.  Ultimately,  the  desire  is  to  improve  the  detectability  of  the 
target  echo.  WiUi  this  in  mind,  little  loss  is  tolerated  in  the  target  echo  level,  while  as  much 
reverberation  is  suppressed  as  possible.  In  this  manner,  the  target’s  signal-to-background  ratio  is 
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enhanced  and  the  echo  becomes  more  detectable.  Now,  the  theoretical  requirements  of  target 
motion  (Doppler)  for  echo  preservation  in  the  subtraction  process  will  be  discussed. 

Consider  the  target’s  motion  relative  to  the  sonar  system,  as  shown  in  figure  3.  The  target  is 
assumed  traveling  with  a  constant  velocity,  vy,  and  constant  heading  over  the  interval  of  time  that 
reverberation  is  being  processed  for  suppression.  Assume  that  the  source  and  the  receiver  are 
fixed  in  location,  with  no  motion.  The  velocity  of  the  target  with  respect  to  a  bistatic  geometry,  v, 
is  givenl^l  as 

V  =  V j.  •  (cos  0„  +  cos  0,;.) 

=  2vj  '  cos0,^  •  cos^ 

where  0;5,  0„.,  Qtb,  ^nd  P  are  various  angles  associated  with  the  geometry  as  shown  in  figure  3.  p 
is  defined  as  the  opening  angle  between  the  source,  target,  and  receiver.  0,*  is  the  angle  from  the 
target  heading  to  the  bisector  of  P  and  is  known  as  the  target  bistatic  aspect  angle.  The  angles  0fj 
and  Qtr  are  defined  as  the  angles  from  the  target  heading  to  the  source  and  receiver,  respectively. 
Notice  that  if  p  is  equal  to  zero  and  4/  =  the  system  is  monostatic,  with  the  source  and  re¬ 
ceiver  colocated.  In  addition,  if  Qft,  =  0,  the  target  is  moving  radially,  with  an  closing  velocity 
component  given  by  2vy. 

Assuming  the  use  of  certain  broadband  signals  processed  with  matched  filtering,  the  Doppler 
of  the  received  target  echo  is  not  observable  as  a  shift  in  frequency  but  as  a  change  in  range  from 
one  ping  to  the  next.  This  change  in  bistatic  range  corresponds  to  the  location  of  the  target  at  dif¬ 
ferent  equitime  ellipses  from  ping  to  ping.  After  alignment  of  the  w'*'  and  n'**  pings  for  subtrac¬ 
tion,  the  time  shift  between  echoes  is  given  by  the  difference  in  arrival  time  relative  to  their 
respective  direct-blast  receptions.  The  time  shift  between  echoes,  At^,  is  given  as 

Ate  =  -  te„  =  — ,  where  AT  =  (m  -  n)  •  (32) 

Assuming  that  the  received  reverberation  is  caused  by  scatterers  that  are  stationary,  and  that 
it  is  coherent  from  ping  to  ping,  it  may  be  suppressed  according  to  the  appropriate  RPR  for  the 
method  used.  To  do  this,  the  parameters  of  the  adaptive  filtering  algorithm  must  be  so  tuned  that 
the  algorithm  is  updating  quickly  enough  to  track  and  cancel  the  distributed  reverberation  back¬ 
ground,  but  not  so  quickly  that  it  cancels  uncorrelated  random  noise.  There  may  or  may  not  be 
subtraction  of  the  input  target  echo,  depending  on  the  amount  of  time  shift  between  the  echoes 
from  ping  to  ping  after  alignment.  At  the  output  of  the  replica  correlator,  in  the  absence  of  noise, 
the  target  echo  is  approximated  by  the  autocorrelation  function  of  a  band-limited,  flat  spectral 
signal.  The  frequency  modulated  chirp  waveforms  have  these  frequency  characteristics.  A 
real-valued  representation  of  a  target  echo  for  these  types  of  waveforms  is  givenf*^!  by 

e(t)  »  cos(2jifct)  =  sinc(;r5r)  ■  cos(2jtfct)  (33) 

where  fc  and  B  are  the  transmitted  signal  center  frequency  and  bandwidth,  respectively  The  echo 
has  a  carrier  frequency  which  is  modulated  by  a  sine  function.  Figure  15  shows  an  example  of  a 
noise-free  echo  signal  at  the  output  of  the  correlator  with  the  ratio  of  center  frequency  to  band¬ 
width  equal  to  3.33.  The  maximum  peak  of  the  target  echo  that  stands  up  above  the  background 
level  provides  the  target’s  detectability.  Therefore,  the  observable  echo  level  at  the  output,  is 
the  maximum  of  the  difference  between  two  identical  echoes,  one  shifted  relative  to  the  other,  as 


Bout  ~  MAX{6|||(f)  —  Cb(/ 

=  MAX'^sinc(;rfir)  •  cos(2nr/c) 

-sincj^;rJ5^r  -  *  cos^27r/c|/  - 

The  input  echo  level,  einput,  is  the  maximum  of  the  single  target  echo  in  the  input  time  series 
which  in  this  case  is  equal  to  one.  The  echo  loss  (EL),  which  is  the  loss  in  echo  power  can  then 
be  given  by  the  ratio  of  input  to  output  levels  as 

“  =  ‘•>'<>810(1^1')  =  (35) 


Figures  16  through  19  show  four  plots  of  EL,  shown  as  a  function  of  echo  time  shift,  A4,  for 
waveform  center-frequency-to-bandwidth  ratios  of^/B  =  1.67, 3.33, 6.67,  and  13.33.  Observe 
that  the  output  echo  power  oscillates  at  the  carrier  frequency,)^,  but  is  bounded  by  an  upper  and 
a  lower  envelope.  Depending  on  the  echo’s  time  shift,  the  output  echo  power  can  undergo  signif¬ 
icant  enhancement  or  cancellation,  particularly  at  small  shifts.  For  values  of  echo  time  shift  less 
than  l/B,  severe  cancellation  and  significant  enhancements  are  realizable.  The  echo  levels  in  this 
region  are  highly  unpredictable  and  therefore  deemed  unreliable.  Clearly,  as  the  echo  time  shift 
increases  beyond  l/B,  there  are  still  enhancement  and  cancellation  peaks,  yet  the  excursions 
from  the  input  echo  level  (0  dB)  are  reduced,  as  they  attenuate  with  increasing  shift.  If,  the  target 
motion  is  such  that  it  yields  an  echo  shift  that  is  in  this  region,  the  target  echo  will  be  reliably 
preserved. 
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EL  (dB)  EL  (dB) 


ECHO  TIME  SHIFT  (s) 

Figure  16.  Echo  loss  vs.  echo  time  shift  {fclB  -  1.67). 


ECHO  TIME  SHIFT  (8) 

Figure  17.  Echo  loss  vs.  echo  time  shift  {fclB  -  3.33). 
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Figure  18.  Echo  loss  vs.  echo  time  shift  (^/B  =  6.67). 


In  order  to  validate  that  these  curves  accurately  reflect  the  behavior  of  the  echo  loss,  target 
echoes  were  synthesized  and  a  simulation  was  performed.  For  this  simulation,  two  time  series 
(pings)  were  produced  with  an  identical  random  noise  signal,  in  order  to  model  a  perfectly 
coherent  reverberation  background.  A  hyperbolic-frequency  modulated  (HFM)  waveform  was 
then  added  to  each  of  the  two  time  series,  as  a  simulation  of  the  echo.  The  echo  waveforms  had  a 


25 


selectable  Doppler  shift  applied  to  them  before  they  were  summed  with  the  reverberation.  This 
Doppler  shift  was  determined  by  specifying  both  the  bistatic  velocity  and  the  ping  repetition 
interval,  as  in  Eq.  32.  The  simulated  echoes  were  given  such  large  negative  signal-to-background 
levels  that  they  were  deeply  buried  in  and  obscured  by  the  simulated  reverberation.  The  two  time 
series  were  then  matched-filtered  with  a  replica  of  the  HFM  waveform.  Even  with  the  resultant 
gain  of  the  matched  filter,  the  target  echoes  were  still  obscured  by  the  reverberation. 

A  subtraction  algorithm  was  then  applied  to  the  two  time  series  and  the  output  echo  levels 
were  measured  from  the  suppression  output.  In  this  manner,  the  echo  loss  was  measured  for  a 
particular  Doppler  assumption.  This  procedure  was  then  repeated  many  times,  with  a  large  set  of 
different  Doppler  assumptions,  in  order  to  trace  EL  as  a  function  of  target  Doppler.  Figure  20 
shows  the  echo  loss  measurements  for  a  set  of  400  subtractions  with  different  Doppler  assump¬ 
tions.  For  this  simulation,  the  ping  repetition  interval  was  set  at  60  seconds,  while  the  bistatic 
velocity  was  varied  over  the  range  from  0  to  6  knots,  so  that  the  echo  time  shift  is  only  a  function 
of  the  target  velocity.  The  embedded  target  echo  waveform  had  a  center-frequency-to-bandwidth 
ratio  of  3.33.  Here,  a  straight  coherent  (nonadaptive)  subtraction  was  performed.  It  is  clear  that 
the  simulation  of  echo  loss  agrees  very  closely  with  what  is  predicted  by  the  theory  (Eq.  34  and 
35).  Next,  the  simulation  was  performed  using  the  adaptive  RLS  filter  for  the  subtraction.  Figure 
21  shows  the  result  of  this  echo  loss  simulation  over  the  same  range  of  velocities  as  before.  Here 
again,  that  the  simulated  echo  loss  is  behaving  as  predicted  by  theory.  Also,  that  the  measured 
echo  loss  of  the  RLS  simulation  is  virtually  identical  to  the  simulation  of  the  nonadaptive  sub¬ 
traction.  The  results  of  this  simulation  indicate  that  the  theoretical  echo  loss  given  in  Eq.  34  and 
35  is  an  accurate  predictor  of  the  type  of  echo  losses  that  will  be  suffered  with  real  target  echo 
data. 


Figure  20.  Echo  loss,  theory  vs.  simulation  (nonadaptive). 
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Now  consider  only  the  two  bounding  envelopes,  which  determine  the  regions  of  echo  shift 
for  reliable  target  preservation.  The  envelopes  themselves  are  solely  a  function  of  the  transmitted 
pulse  bandwidth.  Figures  22  through  25  show  both  envelopes  for  values  of  5  =  10, 20, 30,  and 
40  Hz,  respectively.  Notice  that  as  the  bandwidth  increases,  the  envelopes  attenuate  much  earlier. 
Therefore,  increases  in  bandwidth  allow  for  searches  of  targets  with  smaller  Doppler  shifts,  as 
the  reliable  target  preservation  region  extends  down  to  slower  velocities. 

Now  consider  a  speciflc  example  of  EL  as  shown  in  figure  26.  The  time  shift  is  composed  of 
the  velocity  component  and  the  time  between  pings.  Here  the  bandwidth  is  equal  to  30  Hz,  and 
the  horizontal  axis  is  the  radial  velocity  component,  v,  of  the  target  in  knots.  Several  sets  of 
envelopes  are  shown  overlaid,  for  different  values  of  time  between  pings,  AT.  As  AT  is  length¬ 
ened,  the  reliable  region  for  echo  preservation  is  extended  down  to  slower  Doppler  shifts. 

Figure  27  shows  a  similar  plot  for  different  target  velocities.  Here,  for  a  fixed  time  between 
pings,  faster  moving  targets  are  the  most  reliably  detected.  These  plots  aid  in  the  determination 
of  the  range  of  target  Doppler  shifts  that  will  be  passed  through  the  subtraction  with  minimum 
loss.  In  addition,  they  facilitate  the  determination  of  optimum  values  of  A7  for  detection  in 
certain  target-motion  scenarios.  While  it  is  desirable  to  extend  AT  to  obtain  reverberation 
suppression  for  slow  targets,  there  may  be  limits  to  how  much  this  can  be  done.  It  is  vital  that  the 
interping  coherence  remain  high  over  the  time  period  considered.  For  some  acoustic  environ¬ 
ments,  the  reverberation  field  may  only  be  stationary  over  some  finite  length  interval  of  time.  In 
such  an  event,  the  interping  coherence  may  decreases  in  level  as  AT  is  increased  beyond  this 
time  interval,  lowering  the  obtainable  RPR. 
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Figure  23.  Echo  loss  envelopes  (B  =  20). 
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ECHO  TIME  SHIFT  (s) 

Figure  24,  Echo  loss  envelopes  {B  =  30). 


ECHO  TIME  SHIFT  (•) 

Figure  25.  Echo  loss  envelopes  (B  =  40). 
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Figure  26.  Echo  loss  envelopes  for  varying  repetition  rates. 


PING  REPETmON  II^RVAL  (min) 

Figure  27.  Echo  loss  envelopes  for  varying  target  Dopplers. 


30 


GAIN  AGAINST  REVERBERATION 


Now  consider  the  total  gain  against  reverberation  (GAR)  as  a  combination  of  the  RPR  and  all 
processing  losses.  This  is  the  total  theoretical  gain  in  detectability  for  a  target  echo  in  a 
distributed  reverberation  background,  processed  using  interping  suppression  techniques.  Having 
achieved  this  gain,  the  data  can  then  be  put  through  the  detector,  with  a  better  probability  of 
detection.  The  gain  against  reverberation  for  Case  1  is  given  by  the  reverberation  power  reduc¬ 
tion  for  Case  1,  and  the  target  echo  loss  as 

GARi  =  RPR^  -  EL  (36) 

The  gain  against  reverberation  for  Case  2  is  similarly  obtained,  except  that  it  includes  an 
additional  loss  term,  RL.  It  is  given  as 

GARj  =  RPRi  -  RL  -  EL  (37) 

RL  is  the  reference  loss  due  to  the  input  target  echo  being  included  in  the  reference  signal,  and 
therefore  only  applies  to  Case  2.  Because  the  target  echo  is  in  both  the  input  and  the  reference,  it 
is  partially  subtracted  out  in  the  reverberation  processing.  So  for  this  case,  the  output  target  echo, 
Eout>  is  equal  to  the  input  echo,  einput^  minus  the  fraction  of  the  input  echo  that  went  into  the 
averaged  reference,  so  that 


For  detection  processing,  the  signal’s  power  is  used.  The  square  of  the  output  echo  gives  the 
echo  power  as 


P 


(39) 


With  this  result,  the  loss  due  to  the  reference  can  be  determined  as  the  ratio  of  echo  input  to 
output  power,  in  decibels,  as 


(40) 


Using  Eq.  25  and  40,  the  first  two  terms  of  Eq.  37  are  simpliEed  as 

N 


RPR2  -  RL  =  10  log 


N  .  (N-  1\^ 

-  i)(i  - r)  \  N  ) 


=  101ogio[^(l  -  r)]  =  RPR, 

This  result  shows  that  the  overall  gain  associated  with  Case  1  and  Case  2  is  really  equivalent,  for 
N  total  input  pings.  Thus,  while  Case  2  suppresses  more  reverberation  than  Case  1,  it  also  suffers 
an  additional  loss  in  the  target  echo.  The  increased  reverberation  suppression  exactly  equals  the 
additional  target  loss,  so  they  offset  one  another.  Therefore,  the  final  result  obtained  by  using 
either  the  method  in  Case  1  or  Case  2  is  a  gain  in  detectability  given  by 

GAR  =  RPRy  -  EL  (42) 
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VALIDATION  WITH  DATA 


Now  actual  reverberation  data  collected  by  a  bistatic  receiver  off  the  western  coast  of  the 
United  States  are  used  to  validate  the  interping  coherence  measurement  as  a  predictor  of  suppres¬ 
sion  performance.  The  target  loss  effect  is  not  validated  with  the  data,  because  it  does  not  contain 
any  target  echoes.  The  receiver  was  fixed  on  the  bottom  and  located  at  the  base  of  the  rising  con¬ 
tinental  slope.  The  source  was  to  the  east,  transmitting  a  variety  of  broadband  signals  in  several 
different  frequency  bands.  The  acoustic  transmissions,  as  well  as  backscattered  reverberation 
from  the  rising  continental  slope,  were  received  by  the  bistatic  array.  The  received  signals  were 
processed  using  an  active  sonar  signal  processing  software  package!^!,  in  a  fashion  similar  to 
that  described  in  the  Active  Sonar  Signal  Processing  section.  The  received  data,  the  output  of 
coherent  reverberation  subtraction,  measurements  of  interping  coherence,  and  coherence  of 
reverberation  will  now  be  shown. 

Figure  28  shows  five  sequentially  received  pings,  for  a  single  beam  looking  at  the  continental 
shelf.  The  data  have  also  been  matched  filtered  and  smoothed,  with  the  magnitude  of  the  signal 
power  plotted.  The  received  data  stream  has  been  segmented  from  five  contiguous  blocks  of 
duration  Tpri  =  90  seconds  and  stacked  vertically.  Stacked  in  this  manner,  the  pings  line  up,  with 
the  received  direct  blast  located  at  around  8  seconds.  Although  the  received  signals  are  not  cali¬ 
brated  to  sound  pressure  levels,  the  relative  levels  show  the  received  signal  power  fluctuations. 
Following  the  direct  blast  are  several  areas  of  significant  reverberation  power,  most  noticeably  at 
around  70-85  seconds.  This  reverberation  is  backscatter  from  the  rising  continental  slope.  Each 
of  the  five  received  pings  shows  very  similar  power  characteristics  through  time.  Also  notice 
how  the  ping  repetition  interval  was  chosen  to  be  sufficiently  large  for  the  reverberation  to  die 
away  before  the  next  received  ping  comes  in. 

The  interping  coherence  function  measurement  between  ping  1  and  ping  2  is  shown  in  figure 
29c.  From  these  data  it  can  be  seen  that  the  interping  coherence  is  quite  high  near  the  direct  blast 
and  the  shelf  reverberation  regions.  In  the  other  regions,  it  is  lower,  because  the  reverberation 
power  is  lower  relative  to  the  underlying  noise  background.  Also  shown  is  ip,  the  delay  lag  of 
the  peak  correlation.  It  can  be  seen  that  where  the  reverberation  strongly  dominates  the  noise,  the 
delay  values  are  constant,  in  this  case  at  -0.02  second.  Where  the  reverberation  is  not  as  strong 
the  coherence  values  drop,  and  the  delay  values  of  the  peak  are  more  randomly  distributed,  as 
expected.  From  the  measurement  of  interping  coherence,  the  largest  gains  against  reverberation 
are  expected  to  be  near  the  direct  blast  and  the  shelf  regions  of  the  ping. 

The  interping  coherence  functions  between  ping  1  and  all  the  others  in  the  set  are  shown 
overlaid  in  figure  30.  The  same  basic  trend  is  clearly  shown  by  all  the  plotted  coherence  func¬ 
tions.  In  particular,  the  direct  blast  is  shown  to  be  highly  coherent  at  the  same  level  between  all 
pairs  of  pings.  In  addition,  the  shelf  reverberation  region,  as  well  as  other  distinguishable  rever¬ 
berators,  shows  the  same  structure.  This  implies  that  the  interping  coherence  function  is  rela¬ 
tively  stationary  over  at  least  57^„.  Hence,  for  a  two-ping  subtraction  process,  equivalent  gains 
against  reverberation  can  be  obtained  by  using  a  reference  ping  just  before  or  five  pings  before 
the  input  ping.  This  stationarity  also  allows  consideration  of  a  multiping  subtraction  up  to  at  least 
the  five  pings  shown. 
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Figure  29.  Interping  coherence  measurement,  ping  1  and  ping  2. 
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Figure  30.  Interping  coherence  functions. 
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Figure  31.  Theoretical  reverberation  power  reduction. 


From  the  interping  coherence  function  measured  between  pings  1  and  2,  the  theoretical  RPR 
was  calculated  by  using  Eq.  20  and  25.  The  results  are  shown  in  figure  31  for  two  cases.  Case  1 
is  for  a  suppression  process  with  two  pings,  where  the  input  is  not  a  part  of  the  reference.  Case  2 
is  for  a  process  with  five  pings,  with  the  input  included  in  the  reference  estimate.  The  theoretical 
RPRs  near  the  direct  blast  and  the  shelf  regions  show  that  significant  reductions  are  possible. 
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with  modest  reductions  elsewhere.  The  multiping  process  shows  larger  suppression,  as  expected. 
The  values  of  RPR  are  shown  clipped  at  zero,  because  in  reality  the  adaptive  algorithms  are 
constrained  so  that  they  do  not  yield  higher  power  in  the  output  than  the  power  in  the  input. 
These  reduction  levels  predicted  by  the  interping  coherence  measure  are  now  verified  by  obtain¬ 
ing  the  actual  reduction  levels,  using  adaptive  suppression  algorithms. 

The  data  were  processed  using  two  different  coherent  reverberation-suppression  algorithms. 
The  first  is  a  two-ping  coherent  subtraction  process  with  the  input  ping  not  included  in  the  refer¬ 
ence  (Case  1,  =  2).  The  input  and  reference  are  pings  1  and  2,  respectively,  and  the  RLS  adap¬ 

tive  algorithm  was  used  with  appropriate  choices  of  parameters.  In  the  second  process,  the  input 
(ping  1)  and  the  other  four  pings  were  input  into  an  SVD  algorithm  (Case  1,N  =  5).  Figure  32 
shows  the  smoothed  power  levels  of  the  input  ping,  along  with  the  suppression  outputs  for  both 
methods,  overlaid.  It  can  be  seen  that  the  reverberation  power  levels  are  reduced  most  noticeably 
in  the  area  of  the  direct  blast  and  the  shelf  reverberation.  The  multiping  method  shows  more 
reverberation  suppression,  as  expected. 

To  quantify  the  amount  of  reverberation  reduction  that  these  processes  have  actually 
achieved,  the  power  levels  of  the  suppression  outputs  were  subtracted  from  those  of  the  input. 
This  gives  the  actual  RPR  obtained  by  using  the  algorithms,  as  shown  in  figure  33.  The  actual 
achieved  RPR  can  then  be  plotted  along  with  the  theoretical  RPR,  which  was  obtained  using  the 
measured  interping  coherence.  This  is  shown  in  figures  34  and  35  for  cases  1  and  2,  respectively. 
For  both  cases,  the  theoretical  RPR  closely  matches  the  actual  obtained  RPR.  In  this  manner, 
using  real  reverberation  data,  the  interping  coherence  function  has  been  validated  as  an  accurate 
predictor  of  reverberation  reduction. 

Finally,  the  amount  of  coherence  in  the  reverberation  component  was  estimated  by  removing 
from  the  interping  coherence  the  degradation  caused  by  the  noise.  Figure  36a  shows  the  received 
reverberation  plus  noise  level  for  ping  1.  Also  shown  is  the  estimated  ambient  noise  level,  which 
is  assumed  to  be  constant  over  the  duration  of  the  processing.  From  this,  estimates  of  the  RNNR 
and  Tmax  were  obtained  by  using  Eq.  26  and  27,  respectively.  Tmax  is  shown  in  figure  36b,  along 
with  the  measured  interping  coherence  function  between  pings  1  and  2.  Note  that  there  are  times 
when  virtually  all  the  possible  coherence  is  being  obtained,  as  well  as  times  when  less  coherence 
is  achieved.  Most  noticeably,  the  direct  blast  is  as  coherent  as  is  possible,  while  the  shelf  rever¬ 
beration  is  not,  even  though  their  received  power  levels  are  comparable.  For  the  reverberation- 
limited  case,  only  coherence  values  above  50%  are  considered.  From  these  coherence  values  the 
coherence  of  the  reverberation,  Freverb^  is  obtained  by  using  Eq.  30.  This  is  shown  in  figure  36c. 
Notice  that  the  received  direct  b'ast  is  97%  coherent,  while  the  reverberation  from  the  shelf  is 
80-90%  coherent.  Also,  notice  the  reverberators  at  around  53,  60,  and  65  seconds.  Although 
they  have  much  smaller  power  in  the  received  data,  they  are  more  coherent  than  the  shelf  rever¬ 
beration.  Therefore,  this  shows  that  the  interping  coherence  is  determined  not  only  by  the 
RNNR,  but  also  by  the  coherence  of  the  reverberation.  In  this  case,  it  also  shows  that  it  is 
possible  to  experience  reverberation  with  varying  power  levels  that  have  differing  coherence 
values. 
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Figure  33.  Reverberation  power  reduction  achieved. 
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(b)  Maximum  coherence  obtainable  vs.  measured  interping  cohererKe. 
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(c)  Coherence  of  reverberation. 


Figure  36.  Estimation  of  coherence  of  reverberation. 
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SUMMARY 


We  have  described  the  approach  of  applying  adaptive  interference  cancelling  techniques  to 
active  bistatic  sonar  data,  where  reverberation  is  limiting  the  detection  performance.  The 
performance  criteria  for  both  the  target  echo  and  the  reverberation  have  been  described.  The 
theoretically  expected  target  echo  losses  were  validated  with  a  simulation  of  synthetically 
generated  target  echoes.  The  interping  coherence  measure  has  been  shown  to  be  an  effective 
means  of  predicting  the  bounds  of  reverberation  suppression.  This  has  been  validated  by 
applying  the  interping  coherence  measure  on  actual  reverberation  data  and  comparing  the 
predicted  RPR  with  the  achieved  RPR.  This  comparison  showed  very  good  agreement, 
indicating  that  the  adaptive  algorithms  were  effectively  tuned  to  obtain  the  maximum  possible 
performance  against  reverberation. 
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